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Evidence for multiple superconducting gaps in optimally doped BaFei . 8 7Co .i3As 2 from 

infrared spectroscopy 
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We performed combined infrared reflection and ellipsometry measurements of the in-plane optical 
response of single crystals of the pnictide high temperature superconductor BaFe1.87Coo.13 As2 with 
T c = 24.5 K. Our experimental data provide evidence for multiple superconducting energy gaps 
and can be well described in terms of three isotropic gaps with 2A/fcsT c of 3.1, 4.7, and 9.2. The 
obtained low-temperature value of the in-plane magnetic penetration depth is 270 nm. 

PACS numbers: 74.70.-b,78.30.-j, 74.25. Gz 
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The origin of the recently discovered high temperature 
superconductivity (HTSC) in the pnictides is a subject 
of great interest^ In the pnictides, it has been predicted 
that up to five bands are crossing the Fermi-level and thus 
can participate in superconducting (SC) stated There- 
fore they are likely multi-band superconductors where 
the magnitude of the energy gaps and even the sign of the 
order parameters can vary between different bandsi^r— 
Despite of the tremendous research efforts, there is an 
ongoing debate about the basic superconducting proper- 
ties, such as the number of energy gaps and their mag- 
nitudes, which are prerequisites for determining the or- 
der parameter symmetry and ultimately the HTSC pair- 
ing mechanising— Experimental evidence for at least 
two different gaps has already been reported from angle- 
resolved photoemission (ARPES), point contact spec- 
troscopy, nuclear magnetic resonance, and muon spin ro- 
tation (/^SR)i2r— However, the reported gap values ex- 
hibit a large variation from 2A/fcsT c m 1.6 to 10 that 
remains to be understood. Possible factors are: a varia- 
tion between the different compounds and as a function 
of doping and structural changes, a strong dependence on 
the sample (surface) quality, or the presence of even more 
than two gaps which the various experimental techniques 
are probing with different sensitivity. More detailed and 
accurate measurements of the properties of the SC energy 
gap(s) are thus required. 

Infrared (IR) spectroscopy is a powerful technique to 
investigate the electronic gap structure of superconduc- 
tors. While it does not provide fc-space resolved infor- 
mation such as ARPES, its large probe depth ensures 
the bulk nature of the measured quantities and its high 
energy resolution and powerful sum rules enable a re- 
liable determination of important physical parameters, 
such as the gap magnitude and the plasma frequency of 
the SC condensate. In the parent compound of the Ba 
122 phase, (Ba,Sr)Fe2As2, the IR data revealed the spin 
density wave gap(s) and related phonon anomalies i^Sr— 
For the SC single crystals the first report of the gap fea- 
ture has been published for hole doped Bai_ :E K. c Fe 2 As 2 
with T c — 37 Ki^i In addition to a steep absorption edge 
at 150 cm -1 that is characteristic of a nearly isotropic 



(nodeless) gap with 2A//csT c « 6, they also found evi- 
dence for a gap with 2A/fcsT c « 8. A similar upper value 
of 2A/ksT c R3 8 has been obtained from an ellipsome- 
try study on polycrystalline (Nd,Sm)FeAsOi_ a; F 2 ;. 1? Re- 
cent studies on electron doped BaFe 2 - a; Co : , : As2 single 
crystal s 19 " 21 also reported pronounced gap edges around 
30 - 50 cm . These reports still vary concerning the 
number and the magnitude of the SC energy gaps and of 
the residual low-frequency conductivity. 

In the following we present a detailed IR spectroscopy 
study of electron doped BaFei.87Coo.i3As2 single crystals 
with T c = 24. 5K. Our spectra reveal several SC-induced 
features which yield further insight into the multigap na- 
ture of SC. In particular, we show that a BCS-type model 
with three isotropic energy gaps with 2A/ksT c ps 3.1, 
4.7, and 9.2 describes the data rather well. 

The BaFei.87Coo.i3As2 single crystals were grown from 
self-flux in glassy carbon crucibles and their chemical 
composition was determined by energy dispersive x-ray 
spectroscopy as described in Ref. 9. A bulk SC tran- 
sition of T c = 24.5 K was confirmed by transport, dc 
magnetization, and /iSR measurements. 

The optical measurements were performed on freshly 
cleaved pieces from a single growth batch. The temper- 
ature (T) dependent near normal incidence reflectivity 
spectra, R(co), for 35-5000 cm -1 were measured with a 
Brukcr 113v FT-IR spectrometer utilizing the in situ gold 
evaporation technique^ Additional ellipsometry mea- 
surements were performed at 350-8000 cm -1 with a home 
built rotating-analyzer ellipsometer attached to a Brukcr 
113v (Ref. 23) and at 6000 - 52000 cm" 1 with a commer- 
cial Woollam VASE ellipsometer. All measurements were 
performed at least twice to ensure their reproducibility. 
The ellipsometry data were converted to normal inci- 
dence reflectivity data. The complex optical constants 
were then calculated by Kramers-Kronig transformation 
with proper extrapolations such that the directly mea- 
sured ellipsometric data were reproduced. 

Figure 1(a) shows the reflectivity, R(u>), in the far- 
infrared (FIR) region at selected temperatures. The real 
part of the optical conductivity, ai(u)), is displayed in 
Fig. 1(b). The overall shape and T- dependence of 
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FIG. 2: Self-consist determination of the SC plasma fre- 
quency, u Pi sc- The measured ef K (symbols) is shown to 
compare well with the expected spectrum (solid line) that 
has been derived from the measured 
missing inductive response due to Aai 
set) to obtain e reg (dotted line) and adding the contribution 
e s (dashed line) of the SC delta function. 



£i° K by subtracting the 
af K -a? K (see the in- 



FIG. 1: Temperature dependence of the FIR reflectivity, R(cj) 
(a) and optical conductivity ai (lj) (b) . Solid circles show a dc 
as reported for optimally doped BaFei_ x Co x As2^. 



the spectra are similar to previous reports i ' 19 ~ 21 ' i 25 
In particular, the spectra contain a pronounced Drudc- 
peak at low frequency due to free carriers, which nar- 
rows significantly as T decreases in the normal state. 
The extrapolated dc values, <r dc , compare well with pub- 
lished transport data (solid circles).— The electronic re- 
sponse contains at least two more broad bands in the 
mid-infrared (MIR) range that likely arise from inter- 
band transitions . 14 ' 25 In addition, our data contain two 
sharp features at 95 and 261 cm -1 which correspond to 
IR-active phonons. In the SC state at T < T c , the elec- 
tronic response undergoes some characteristic changes 
below about 300 cm -1 which provide detailed informa- 
tion about the SC energy gap(s) and the SC condensate. 

First we discuss the determination of the plasma fre- 
quency of the SC condensate, u> P: sCi which is detailed 
in Fig. 2. The inset shows the SC-induced decrease in 
the regular part of the conductivity, A<Ji(uj)=af OK (oj)- 
from which we obtain the missing spectral 
weig ht, AS* = f™Aai(co)duj. Its magnitude falls off 
steeply with increasing frequency and essentially vanishes 
above 500 cm -1 (see the inset of Fig. 2). According 
to the so-called Ferell-Glover-Tinkham sum rule, AS is 
redistributed to a delta function at the origin which ac- 
counts for response of the SC condensate^ The analysis, 
using a low-frequency extrapolation of the 30 K spec- 
trum of er dc = 10000 il^cm -1 as reported in Ref. 26, 
yields AS « 9.1 x 10 5 f2 -1 cm~ 2 corresponding to 0J Pt sc ~ 
5900 cm -1 and a value of the magnetic penetration depth 



of X a b = - c sc ~ 270 nm. This agrees reasonably well 
with Xai, R3 230 nm as obtained from /iSR measurements 
on optimally doped BaFei_ 2: Coa;As2.— ^ The self consis- 
tency of this analysis has been checked by inspecting the 
corresponding changes in the real part of the dielectric 
function, ei(ui). Figure 2 shows the comparison between 
the measured ef K (w) (symbols) and the corresponding 
spectrum (solid line) that is expected based on the SC- 
induced spectral weight redistribution and, in particu- 
lar, the low frequency extrapolation that was used to 
obtain uj p .sc- Specifically, from £ 30K (cj) we subtracted 
the missing contribution of the regular response due to 
Aai(oj)=af° K (oj)-af K (uj) to obtain e reg (dotted line) and 
added the contribution e s (dashed line) due to the SC 
delta function. The excellent agreement between the 
symbols and the solid line confirms that we performed 
a reasonable low- frequency extrapolation of o\ . 

Next we show that valuable information about the 
SC energy gaps can be obtained from the SC-induced 
changes in <J\(iS) and R(lu). First, we note that the pro- 
nounced gap edge around 50 cm -1 is the hallmark of an 
(almost) isotropic energy gap whose magnitude does not 
far exceed the normal state scattering rate, 7. In that 
case cri(w) exhibits a pronounced feature at 2 A where it 
increases sharply from zero at w < 2 A to a maximum 
around 4A, before it eventually merges with the nor- 
mal state data at higher frequency . 30 ' 31 Correspondingly, 
R{uj) remains unity at lj < 2 A and suddenly starts to 
decrease at tu > 2A yielding a characteristic maximum 
in the reflectivity ratio R(T < T C )/R(T > T C )M In a 
multigap superconductor, the corresponding features due 
to the larger gaps are superimposed and thus more diffi- 
cult to observe. Nevertheless, as shown in the following, 
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FIG. 3: Simulation with a BCS-type model of multiple 
isotropic gaps. Comparison of the experimental spectra of the 
(a) conductivity and (b)-(c) reflectivity ratio R(T)/R(30K) 
(symbols) with the simulated spectra (lines) as obtained by 
first fitting the normal state response with a Drude-Lorentz- 
model and then imposing one, two, or three independent 
isotropic BCS-type energy gaps (parameters are given in Ta- 
ble I). Arrows mark clear features due to the gap edges at 2A. 
The inset shows the corresponding experimental (symbol) and 
simulated (red line) spectra of the dielectric function. 



the gap magnitudes can still be obtained based on an 
analysis with multiple independent isotropic BCS-type 
superfluids^i 

At first we fitted the normal state spectra of ai° K (oj) 
and ef 0K (w) with a Drude-Lorentz model. Since theory 
predicts that up to five bands are crossing the Fermi- 
level, the low frequency response could be composed of 
as many Drude bands. Nevertheless, we found that two 
Drude oscillators (with parameters as given in Table I) 
are sufficient for a reasonable description of the normal 
state data. Likely these represent the electronlike and 
holclike bands, respectively. Besides, we had to include 
two broad Lorentz oscillators centered at about 1100 and 
5400 cm -1 to account for the electronic MIR bands plus 
two narrow ones at 95 and 261 cm -1 for the IR active 
phonons. Their contribution is shown by the pink line 



in Fig. 3(a) and assumed not to change below T c . Note 
that unlike in Rcf. 19, our data yield no evidence for 
a sizeable electronic mode below about 500 cm -1 other 
than Drude terms^S 

First we demonstrate that the model of a single BCS- 
type gap fails to describe our data. In this case, the sharp 
edge in a\{u), as marked by the black arrow in Fig. 3(a) 
defines the magnitude of the SC energy gap of 2A m 53 
cm -1 . Nevertheless, if we impose this gap value to the 
Drude-bands as obtained at 30 K, the model provides a 
poor description of the data since it largely overestimates 
the magnitude of a\(u) > 2A), in particular, around the 
maximum near 4 A. 

Next we show that a good agreement with the exper- 
imental data can be obtained with a three gap model. 
The signatures of the two larger gaps are considerably 
weaker than the ones of the smallest gap, but they can 
still be identified in the reflectivity ratio, i?(5K)/i?(30K), 
in Fig. 3(b) which eliminates any additional structures 
due to extrinsic features that may be introduced dur- 
ing the thermal cycle in between the gold correction 
procedure. While i?(60K)/i?(30K) shows no anomaly, 
i?(5K)/i?(30K) reveals two noticeable kinks near 80 and 
157 cm -1 as marked by the black arrows. These kinks 
are well separated from the narrow features due to the 
phonon modes at 95 and 261 cm -1 . Therefore they are 
supposed to be SC-induced ones, which should corre- 
spond to the edges of the larger gaps. Concerning the 
assignment of these gaps to the Drude bands (as obtained 
from the 30 K spectrum), the best fits were obtained if 
the gap structures at 53 and 157 cm -1 (80 cm -1 ) were 
associated with the narrow (broad) Drude-band as de- 
tailed in Table I. The quality of the fits was noticeably 
reduced if we changed the assignment of the gaps at 53 
and 80 cm . Even for the largest gap, that is associated 
with a rather small amount of spectral weight, the com- 
parison with the experiment become significantly worse 
if we assigned it to the broad Drude band. The best fit 
with this three gap model is shown by the red lines in 
Fig. 3. It reproduces all the characteristic features of 
the experimental data, including the long tail up to 150 
cm" 1 in i?(5K)/i?(30K)322 

Next we show that the best fits with a correspond- 
ing two-gap model are significantly worse. Figure 3 (a) 
shows the results for three representative assignments (A 
to C) of the gaps to the Drude-bands. It highlights that 
the steep rise of ai(ui) just above 53 cm -1 can only be 
reproduced if the smaller gap at 2 A = 53 cm -1 is as- 
signed to the narrow Drude-band as in models B and C. 
Nevertheless, model B (C) severely underestimates (over- 
estimates) the magnitude of <7i(w) at higher frequency. 
Furthermore, the two gap models fail to reproduce the 
pronounced high frequency tail in i?(5K)/i?(30K) as is 
shown in Fig. 3(c) for model B which yields the most 
reasonable fit with two gaps to <j\(uj). We note that al- 
ternative ways of fitting, for example by describing the 
normal state spectrum with more than two Drude-bands, 
did not improve these two-gap model fits. Also, a strong 
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TABLE I: Parameters for the simulations shown in Fig. 3. 





Drude 1 (7 = 90 cm -1 ) 


Drude 2 (7 = 300 cm' 1 ) 


2 A (cm" 1 ) 


Up (cm z ) 


2 A (cm" 1 ) 


Up (cm *) 


30 K 





4.45x10' 





4.3x10' 


2 gaps-A 


80 


4.45x10' 


53 


4.3x10' 


2 gaps-B 


53 


4.45x10' 


80 


4.3x10' 


2 gaps-C 


53 


4.45x10' 


157 


4.3x10' 


3 gaps 


53 


3.75x10' 


80 


4.3x10' 


157 


7.0x10° 







coupling to a bosonic mode, while it might account for 
one of the kink structures in i?(5K)/i?(30K), would fur- 
ther increase ai(co) above the boson energy and thus 
enhance the discrepancy with the experimental data*22 
Finally, we note that the inclusion of an additional low- 
frequency electronic mode did not help to improve the 
quality of the fits comparable to the one of the three gap 
model— 

Our multigap analysis thus provides evidence for three 
different energy gaps with 2A//csT c w 3.1, 4.7, and 9.2 
in optimally doped BaFe2- :c Co a ;As2. These fall well into 
the range of the reported values of 2A/k B T c w 1.6 - 10 
as obtained with various techniques We note that 
our optical data are compatible with the presence of a 
fourth energy gap of magnitude 2A/ksT c < 3 or with gap 
nodes. As long as this gap or the nodal states involves 
a relatively small amount of spectral weight, we would 
not be able to identify the corresponding signature in 
the spectrum of oi{ui) where the error bars become size- 
able below 53 cm" 1 . Some evidence (though not entirely 
conclusive) for an energy gap with 2A <40 cm -1 or gap 
nodes is contained in the spectrum of i?(5)K/i?(30). As 
shown in Fig. 3(c), it does not exhibit the expected de- 
crease below the lower gap edge at 53 cm -1 , instead it 



remains almost constant down to at least 40 cm -1 . 

Concerning the assignment of these gaps to the differ- 
ent energy bands our optical data provide only limited in- 
formation. This requires additional input concerning the 
momentum dependence for example from ARPES. We 
also note recent theoretical work which considers an or- 
der parameter with extended s-wave symmetry and finds 
that impurity scattering can lead to a substantial modifi- 
cation of the gap structures in the optical conductivity^ 
Clearly, further experimental and theoretical work is re- 
quired to identify the number, magnitude and symmetry 
of the energy gaps in these multi-band superconductors. 

In summary, with combined IR reflection and ellipsom- 
etry measurements we investigated the T-dependent op- 
tical in-plane response of BaFei.87Coo.i3As2 single crys- 
tals with T c = 24.5 K.We showed that the spectra dis- 
play clear SC-induced changes that are characteristic of 
at least three different energy gaps. In particular, we 
showed that they can be well accounted for in terms of a 
BCS-type multigap model with three isotropic gaps with 
2A/ k B T c ss 3.1, 4.7, and 9.2. In addition, we determined 
the SC plasma frequency and the corresponding low-T 
value of the in-plane magnetic penetration depth of 270 
nm. 
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